Oscillating chiral currents in nanotubes: a route to nanoscale magnetic test tubes. 
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With a view to optimising the design of carbon-nanotube (CNT) windmills and to maximising 
the internal magnetic field generated by chiral currents, we present analytical results for the group 
velocity components of an electron flux through chiral carbon nanotubes. Chiral currents are shown 
to exhibit a rich behaviour and can even change sign and oscillate as the energy of the electrons 
is increased. We find that the transverse velocity and associated angular momentum of electrons 
is a maximum for non- metallic CNTs with a chiral angle of 18°. Such CNTs are therefore the 
optimal choice for CNT windmills and also generate the largest internal magnetic field for a given 
longitudinal current. For a longitudinal current of order 10~^ amps, this field can be of order 
10~^Teslas, which is sufficient to produce interesting spintronic effects and a significant contribution 
to the self inductance. 

PACS numbers: 73.63.-b,68.65.-k,71.15.Ap 



Chiral nanotubes and nanowires are of interest for a 
range of properties associated with external magnetic 
fields , their potential for creating nanoscale inductors 
[2] and their suggested role as building blocks in chiral 
nanomotors [3]. Examples studied to date include chiral 
carbon and BC2N nanotubes 0, IH, IgL BN nanotubes 0, 
Fe-filled carbon nanotubes (CNTs) 8| and chiral single- 
wall gold nanotubes [9]. 

Most recently, interest in chiral currents has been 
rekindled by their potential to drive CNT nanometre- 
scale motors [lo|, [ill, 01- Such motors benefit from 
low inter-wall friction [l3| and a high tensile strength 
14 1, which allow one to engineer complex structures 
I5I \m. |l7K including nanoscale bearings [lo|, [l^, ro- 
tors [13, [19 ^oscillators [13, [lH, [H, Q , switches [2^ and 
telescopes 25|. The helical arrangement of the atoms in 
chiral CNTs can be exploited [2^,|23]to produce a Brown- 
ian ratchet effect [sj , rotational and translational motion 
driven by thermal gradients [28, 29] and motion [s^l in- 
duced by circularly polarized light. Recently a new drive 
mechanism for CNT windmills was proposed [3l|, based 
upon the torque generated by a flux of electrons passing 
through a chiral CNT. It was shown that under appropri- 
ate conditions, the dominant contribution to this torque 
is proportional to the flux of angular momentum carried 
by electrons moving in the corresponding infinite chiral 
CNT. This provides a useful guide for the design of CNT 
windmills, since a calculation of the angular momentum 
and associated chiral currents carried by electrons in in- 
finite chiral nanotubes does not require the solution of a 
scattering problem. 

The fact that large chiral currents occur in CNTs is at 
first sight surprising, since early first-principles studies 
[I, 5] suggested that chiral CNTs do not carry a signif- 
icant chiral current. In this paper we study the energy 
and voltage dependence of chiral currents in CNTs and 



show that the small chiral current found in [4] is a con- 
sequence of the metallic nature of the CNT studied. In 
contrast for non- metallic CNTs, we predict much larger 
chiral currents. We find that the energy dependence of 
chiral currents is surprisingly rich. Indeed the transverse 
current components can even change sign and oscillate as 
the energy of the electrons is increased. We demonstrate 
that the presence of large chiral currents produces signifi- 
cant magnetic fields of order 0.1 Teslas within the volume 
of a CNT, thereby providing a novel magnetic test tube, 
which could be used to manipulate the magnetic mo- 
ments of encapsulated magnetic molecules or particles. 
This internal field produces a significant contribution to 
the self inductance of the CNT, which must be added to 
the more usual contribution associated with the external 
magnetic field [32]. 

To define the velocity components of electrons in chi- 
ral CNTs, we follow the notation of [s^, which intro- 
duces the lattice vectors of the corresponding infinite 
2D graphene sheet, defined by ai = ^ and 

a2 = (^-^, — a, where a = \fZac-c and ac-c = 1.44A 
is the carbon-carbon bond length. An (n, m) CNT, where 
(0 < m < n) are integers, is then defined by a transverse 
chiral vector Ch = nai + ma2, which wraps around the 
CNT circumference and a longitudinal translation vec- 
tor T. We are interested in resolving electron velocities 
along axes x and y, which are parallel to the unit vec- 
tors Ch and T respectively. The velocity components 
are given by Hv^ = 9E(k)/9kx and Hvy = 9E(k)/9ky, 
where kx = k.Ch, ky = k.T and E(k) is the energy dis- 
persion relation. In the simplest Slater-Koster scheme, 
E(k) takes the form 

E(k) = 7|1 + exp(-ik.ai) + exp(-ik.a2)|, (1) 

where 7 is the hopping integral. 
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Each mini band possesses a continuous longitudinal 
wave vector ky and is labelled by a quantised value of kx, 
given by 

k5 = ^, (q=l,---Nhex), (2) 

where Nhex is the number of hexagons in a CNT unit 
cell. For a given choice of E and q, equation ([1]) can be 
solved to yield two values of ky. One of these values, 
which we denote k+(q, E), corresponds to a positive lon- 
gitudinal velocity Vy(q, k+(q, E). The other value, which 
we denote k" (q, E), corresponds to a negative longitu- 
dinal velocity Vy(q, k" (q, E). In what follows, we refer 
to these electrons as "right-moving" and "left-moving" 
respectively. We are interested in the transverse ve- 
locities of right-moving electrons, which we denote by 
Vx(q, k+(q, E)). Our aim is to compute the total trans- 
verse velocity Vx"'"^^(E) of all right-moving electrons of 
energy E, which in units of the Fermi velocity vp is 

v(">-)(E) = ^v,(q,k+(q,E))/vF, (3) 

q 

where the sum is over all mini bands with real longitudi- 
nal wave vectors of energy E. 

At low-enough energies, the wave vectors of mini 
band q [namely (k^,k+(q,E)) and (k^, k-(q, E))] can 
be chosen to be close to the K point K. Since E(k) 
is an even function of k, there will be another pair 
of wave vectors in the vicinity of the second K point 
-K, given by (k^', k-(q^ E)) = (-k^, -k+(q, E)) and 
(k^ , k+(q', E)) = — (k^, k" (q, E)), which possess negative 
and positive longitudinal group velocities respectively. 

The contribution to the sum in equation (|3|) from these 
two mini bands is 

Vx(q, E) = Vx(q, k+(q, E)) + Vx(q^ k+(q^ E)) 
= Vx(q, k+(q, E)) + Vx(-q, -k;(q, E)) 

= Vx(q,k+(q,E))-Vx(q,k;(q,E)) (4) 

The last line in this expression is useful, because it 
allows us to focus on the contributions from a single K 
point only. It also demonstrates that a non-zero trans- 
verse velocity arises from trigonal warping, since for a 
perfect Dirac cone, the right hand side of equation (|4]) 
would vanish. This suggests that at low energies, an an- 
alytical expression for Vx(q, E), can be obtained by writ- 
ing the electron wave vector in the form k = K+ry, where 
K = (0, l)47r/(3a) is a vector pointing from the origin to 
a K-point and Taylor expanding E(k) as a power series 
in T^x and r^y, where r]^ = rj.Ch and r]y = ry.T. This 
expansion is of the form 

oo 

= E ^iJ^x^^ (5) 

i,j=0 



where the coefficients cy satisfy cqo = cio = cqi = 
cii and co2 C20 37^/(4a^), C21 

-C03/3 = 277^a^mn(m + n)/16(n^ + m^ +nm)^/^, C12 = 
3V^7^a^(m - n)(2n + m)(2m + n)/16(n^ + m^ + nm)^/^. 
Differentiating this with respect to rj^ and writing 
(k3,k±(q,E)) = K + (7/x,^y ), yields to order [r/^]^] [r/x], 

2Evx(q,E)=ci2[«)'-(7/y-)'] 

+ 7?42C2l[7/+ -7?;] +2C22[«)' - {V^)']} (6) 

To compute r]^ for fixed r]^ and E, we consider two 
cases: The ffist case arises when 77x 7^ 0, in which case 
one obtains 

v.(q,E)/vF = ^^VEVco2 -vi (7) 

The second case corresponds to 77x = 0. In this case, 
equation ([5|) yields to lowest order, E^ = co2^y + co377y, 
which after solving by iteration and combining with equa- 
tion (|6|) yields 

Vx(q,E)/vF = -co3Ci2EVc^f (8) 

where vp = /h = y/3j/2dih is the Fermi velocity. 

Since kx = 27rq/|Ch| and Kx = K.Ch = 27r(n - 
m)/3|Ch|, the value of r]^ for the lowest-energy mini band 
is r]i 27r(-X/3)/|Ch|, where X = 1 for n = m + 1, 
n = m + 4, etc, whereas X = — lforn = m— l,n = m + 2, 
n = m + 5, etc and X = for n = m, n = m±3, etc. Val- 
ues of 77x for higher-energy mini bands are obtained from 
771 by adding or subtracting integer multiples of 27r/|Ch. 
For non- metallic CNTs, where X = ±1, the value of r]^ 
for the second mini band is r]2 = 27r(2X/3)/|Ch|, whereas 
for metallic CNTs, r]2 = ±27r/|Ch|. 

For r]j equation ([7j) yields for the dimensionless 
transverse velocity associated with channel r]j of a (n,m) 
CNT, 

v(-'-)(ryj,E) = Vx(q,E)/vp = v^ef{e - ejj^V^, (9) 
where 

Vj = 3A/6mn(n + m)/ (n^ + m^ + mn)^/^ [sign of t^j] . (10) 

In this expression, e = E/7 and cj is the energy min- 
imum of the jth mini band (in units of 7), given by 
ej = (\/3/2)a|7?j|. 

The above result applies to all low-energy mini bands, 
except the ffist mini band of metallic CNTs, for which 
771 = 0. In this case, equation (|8|) yields 

v(">-)(0,E)=voe2, (11) 

where 

\/3mn((m^ - n^)(2n + m)(2m + n) 

^0 ^ TT^^ 2^ \3 • (12) 

4(n^ + m^ + mn)"^ 
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This quadratic dependence on e means that low-energy 
transverse currents in metahic CNTs is indeed smah, in 
agreement with 0]. In contrast the square-root depen- 
dence arising when 77j 7^ means that transverse currents 
in non- metahic CNTs are predicted to be much larger. 
This behaviour is illustrated in the exact results of Fig.[Tl 
obtained by differentiating equation ([T|) with respect to 
the transverse and longitudinal wave vectors. For each 
mini band q, the red curves of Fig. [T] show the dimen- 
sionless velocities Vx(q, k+(q, E))/vf as a function of E 
for the (8,m) family of CNTs. The black curves show 
the quantity v4"'"^^(E), obtained by adding the values of 
the red curves for each open channel of energy E. As ex- 
pected, for the achiral (8,8) CNT v^^''^\e) = 0, whereas 
for the chiral CNTs y^^''^\e) ^ 0. 
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FIG. 1: Velocity components for the (8,m) family of CNTs 
plotted against the energy E (eV), where m=8,4,5,6 in figs(a) 
to (d) respectively. The red curves show the transverse veloci- 
ties Vx(q, ky (q, E))/vf of right-moving electrons belonging to 
individual channels q. The black curves show the total veloc- 
ity vi"'"^^(E) of Eq. ©. 



As well as predicting the energy dependence of trans- 
verse velocities, equations (p^Oj) and (p!2|) also yield the 
sign of Vx"'"^^(?7j, E). For example, when n and m are 
positive, the sign of vi is equal to the sign of (-X) and 
therefore for the first open mini band of a (8,4) CNT, 
v4^'^^(7^i, E) is negative, whereas for the first mini band 
of a (8,6) CNT, Vx^'^^(77i, E) is positive. Similarly, when 
n > and m > 0, the sign of vq is equal to the sign of 
(m — n). Hence equation ([12]) shows that for the lowest 
mini band of a (8,5) CNT, v4^°''^^(0, E) < 0. For succes- 
sive higher-energy mini bands, r]j alternates in sign and 
therefore equations (|9]) and (p!Q|) reveal that the trans- 
verse velocities of successive higher-energy mini bands 
have a square-root e dependence, with an alternating 
sign. In Eq. (|3]), the label q sums over N(E) right moving 
open channels, where N(E) is a discontinuous function of 
E, which changes by an integer whenever new channels 



open or close. As predicted by equation (|4|), the red 
curves in Fig. (TJ show that right-moving channels open 
or close in pairs and just as a pair of channels open, their 
tangential velocities cancel. Consequently, as shown in 
Fig. [H v4"'"^^(E) is a continuous function of E, with a 
discontinuous first derivative. 

For the purpose of designing a CNT windmill with the 
largest torque, or a CNT with the largest internal mag- 
netic field, it is of interest to compute the maximum in- 
tegrated flux of transverse momentum carried by right- 
moving electrons in an infinite chiral CNT. Since CNTs 
with 771 7^ possess the most favourable energy depen- 
dence for Vx"'"^^/?!, E), we focus on non- metallic CNTs. 



Since 



(n,m) / 



(7^1, E) and v4"'"^^(772, E) have opposite signs. 



the sum Vx"'"^^(e) 



(n,m) / 



(7^i,E)+Vx'''"^^(7^2,E) increases 



monotonically with E for ei < e < 62 and then decreases 
with E for e > 62, passing through zero when 6 — Cmax- 
Hence the maximum integrated flux of transverse veloc- 
ity is proportional to 



/n,m) 



/" 

J ex 



dev("-'")(e), 



(13) 



From equation ([9]), one obtains Cmax = 

1562/14 and 



(n,m) _ /3mn(n + m).[signof (-X)] 



(n2 



m^ 



mn)^/^ 



(14) 



where ^ = 2>/6 7r^ 4arctan^^^^ — arctan ^^^^ ^ 
10.9. 

Equation ([Hj) reveals that Vma^^ ^ as n, m ^ 00, 
which reflects the fact that the angular momentum car- 
ried by an electron wind is a consequence of the finite di- 
ameter of the CNT and the finite difference between suc- 
cessive values of r^j . We also note that the optimum values 
of (n,m), which maximise Vmax^^ are those which possess 
a chiral angle close to 60° (l - f cos"^ v'll/20^ 18°. 

Having analysed the transverse velocity of electrons 
in a chiral nanotubes, we now estimate the magnetic 
field generated by these electrons. In what follows 
we assume that the chiral CNT can be approximated 
by a long solenoid with a constant magnetic field B 
inside the CNT and zero field outside. Since the 
number of right-moving electrons per unit length in 
channel q is (l/7rfi)dE/vy(q, k+(q, E)), and these pass 
around the circumference of the CNT in a time r = 
|Ch|/vx(q, k+(q, E)), the contribution to the tangen- 
tial current per unit length from the qth mini band is 
[(l/7rn)dE/vy(q,k+(q,E))]e/r. The magnetic field B in- 
side a solenoid is /io multiplied by the tangential current 
per unit length. Hence the field due to all right-moving 
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electrons in an energy window eV is 



B = 



2e/io 
hlChl 



Jo 



Vx(q,k+(q,E)) 
Vy(q,k+(q,E)) 



]©(E-Eq)), 



(15) 

where Eq is the lowest energy of the qth mini band. 
Since the current carried by these electrons is I = 
{I/ttH) Jg^^dE^q 6(E — Eq)), the magnetic field can be 
written B = /ioI<^/|Ch|, where 



v^(q,k+(q,E)) 
y(q,k+(q,E)) 



]e(E-Eq)) 



/o^''dEEq©(E-Eq)) 



(16) 



The dimensionless parameter a is the average ratio of 
the tranverse and longitudinal group velocities. At low 
voltages, Eq. (|4j) allows this to be written as a sum over 
channels near the K-point K and using Eq. ([9]), one ob- 
tains 



a = 



/o''deEje(e 



(17) 



For ei < eV/7 < 62, this yields a = viei/\/2, which is of 
order unity for (8,4) or (8,6) CNTs. Consequently, for a 
current of 10~^ amps, |B| ?^ 10~^ Teslas, which is large 
enough to produce significant spintronic effects [35] , such 
as rotating the magnetic moment of a small magnetic is- 
land [36] or a metallocene encapsulated within a CNT 
[37]. By computing the energy stored in this magnetic 
field, we find an associated inductance per unit length 
of L = /io<^^/47r. This internal field could be detected 
through NMR measurements on encapsulated spins and 
could form the basis of scanning MR probe with nanome- 
tre spatial resolution. Finally we note that if the longitu- 
dinal current is driven by an ac voltage, the oscillations 
present in vi"'"^^(E) will lead to the generation of higher 
harmonics, which may provide an alternative probe of 
chiral currents. 
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